Regulator independence of neutron-matter predictions with chiral interactions 
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We calculate the energy per particle in infinite neutron matter perturbatively using chiral N 3 LO 
two-body potentials plus N 2 LO three-body forces. The cutoff dependence of the predictions is 
investigated by employing chiral interactions with different regulators. We find that the inclusion 
of three-nucleon forces, which are consistent with the applied two-nucleon interaction, leads to, 
essentially, regulator-independence of the results. 



PACS numbers: 21.30.Fe,21.65.Cd,21.60.Jz 

I. INTRODUCTION 

A major breakthrough in the last decade has been 
the derivation of nucleon-nucleon (NN) potentials, Vnn, 
based on chiral perturbation theory (ChPT) that are able 
to reproduce accurately the NN data [IH3]. 

The idea to construct realistic two- and three-nucleon 
forces (2NF and 3NF) starting from a chiral Lagrangian 
goes back to the seminal work of Weinberg 4-6 , who 
invoked the concept of an effective field theory (EFT) 
to study the S'-matrix for processes involving arbitrary 
numbers of low-momentum pions and nucleons. In this 
approach, the long-range forces are ruled by the sym- 
metries of low-energy QCD (particularly, spontaneously 
broken chiral symmetry), and the short-range dynamics 
is absorbed into a complete basis of contact terms that 
are proportional to low-energy constants (LECs) fit to 
2N data. 

One great advantage of ChPT is that it generates 
nuclear two- and many-body forces on an equal foot- 
ing [31 El [S]. Most interaction vertices that appear in 
the 3NF and in the four-nucleon force (4NF) also oc- 
cur in the 2NF. The parameters carried by these vertices 
are fixed (along with the LECs of the 2N contact terms) 
in the construction of the chiral 2NF. Consistency then 
requires that for the same vertices the same parameter 
values are used in the 2NF, 3NF, 4NF, .... 

A crucial theme in EFT is regulator independence 
within the range of validity of the theory. In other words, 
the physical observables calculated in the theory must be 
independent both of the choice of the regulator function 
as well as its cutoff scale A. ChPT is a low-momentum 
expansion which is valid only for momenta Q < A x ~ 1 
GeV, where A x denotes the chiral symmetry breaking 
scale. Therefore, NN potentials derived in this frame- 
work are usually multiplied by a regulator function 

f(p',p) = exp[-(p'/A) 2 " - (p/A) 2n ) , (1) 



where typical choices for the cutoff parameter are A ~ 0.5 
GeV. In regards to the physics of the two-nucleon prob- 
lem, it is obvious that the solutions of the Lippmann- 
Schwinger equation, that are related to the two-nucleon 
observables, may depend sensitively on the regulator 
and its cutoff parameter. This unwanted dependence is 
then removed by a renormalization procedure, in which 
the contact terms are re- adjusted to reproduce the two- 
nucleon phase shifts and data. However, it is well-known 
that phase equivalent potentials do not necessarily yield 
identical results in the many-body problem. Thus, one 
may be confronted with cutoff dependence in the many- 
body system [3]. However, in the many-body problem, 
also 3NF, 4NF, ... contribute, which will have impact on 
the final predictions and may either increase or reduce 
the cutoff dependence. 

A convenient theoretical laboratory to investigate this 
issue is infinite nuclear matter and neutron matter. The 
advantage of pure neutron matter is that the contact in- 
teraction, Ve, and the l7r-exchange term, Vd, that ap- 
pear in the N 2 LO three-body force, vanish [TO]. Thus, 
the low-energy constants of Ve and Vd (known as ce 
and Cd), which cannot be constrained by two-body ob- 
servables, are not needed. Consequently, the calculation 
of the ground state energy of infinite neutron matter, 
with chiral 3NFs up to N 2 LO, depends only on parame- 
ters that have been fixed in the two-nucleon system. 

We note that there have been already some attempts 
to study the uncertainties in neutron matter predictions 
using chiral forces, e. g., by Hebeler and Schwenk |10j . 
and Tews et al. |11) . who come up with uncomfortably 
large uncertainties — for reasons to be discussed below. 
It is also worth noting that, aside from the above con- 
siderations, neutron matter, and more generally isospin- 
asymmetric nuclear matter, is currently of great interest 
in the nuclear physics community because of its close con- 
nection with the physics of neutron-rich nuclei and, for 
higher densities, with the structure of neutron stars. 
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It is the purpose of the present paper, to investigate 
how the equation of state of neutron matter, calculated 
using chiral nuclear potentials, depends on the choice of 
the regulator function. More precisely, we employ three 
different chiral potentials, whose cutoff parameters are 
A = 414 [12], 450, and 500 MeV [HE] and calculate, in- 
cluding 3NF effects, the energy per nucleon for neutron 
matter at nuclear densities in the framework of many- 
body perturbation theory. The crucial point of our cal- 
culations is that we use in the 3NF exactly the same 
LECs as in the 2NF. We will show that this consistent 
use of the LECs in the 2NF and 3NF leads to a substan- 
tial reduction of the regulator dependence of the neutron 
matter predictions. 

The paper is organized as follows. In Sec. II, we briefly 
describe the features of the different chiral potentials em- 
ployed and, in Sec. Ill, we give an outline of the calcu- 
lation of the energy per nucleon in neutron matter that 
takes into account 3NF effects. Our results are presented 
in Sec. IV and some concluding remarks and an outlook 
are given in Sec. V. 



II. THE CHIRAL POTENTIALS 

During the past two decades, it has been demonstrated 
that chiral effective field theory (chiral EFT) represents 
a powerful tool to deal with hadronic interactions at low 
energy in a systematic and model-independent way (see 
Refs. [s2 [T3] for recent reviews) . For the construction of 
an EFT, it is crucial to identify a separation of scales. 
In the hadron spectrum, a large gap between the masses 
of the pions and the masses of the vector mesons, like 
p(770) and w(782), can clearly be identified. Thus, it is 
natural to assume that the pion mass sets the soft scale, 
Q r~j m^, and the rho mass the hard scale, A x <~ m p ~ 1 
GeV, also known as the chiral-symmetry breaking scale. 
This is suggestive of considering a low-energy expansion 
arranged in terms of the soft scale over the hard scale, 
(Q / h x ) v , where Q is generic for an external momentum 
(nucleon three-momentum or pion four-momentum) or a 
pion mass. The appropriate degrees of freedom are, ob- 
viously, pions and nucleons, and not quarks and gluons. 
For this EFT to rise above the level of phenomenology, it 
must have a firm link with QCD. The link is established 
by having the EFT observe all relevant symmetries of 
the underlying theory, in particular, the broken chiral 
symmetry of low-energy QCD [3]. The past 15 years 
have seen great progress in applying ChPT to nuclear 
forces. As a result, NN potentials of high precision have 
been constructed, which are based upon ChPT carried 
to N 3 LO. 

Since ChPT is a low-momentum expansion, valid only 
for momenta Q < A x , the potentials are either abruptly 
set to zero for momenta above a certain cutoff A < A x 
("sharp cutoff") or they are multiplied with a smooth 
regulator function, like, e. g., the one of Gaussian shape 
given m Eq. 0. 
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FIG. 1: (Color online) Neutron-proton phase parameters as 
predicted by chiral N 3 LO potentials with different cutoff scale 
A. Solid (red) curve, A = 414 MeV; dashed (blue) curve, A = 
450 MeV; and dotted (black) curve, A = 500 MeV. Partial 
waves with total angular momentum J < 1 are displayed. The 
solid dots and open circles are the results from the Nijmegen 
multi-energy np phase shift analysis [14] and the VPI/GWU 
single-energy np analysis SM99 [15], respectively. 
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FIG. 2: (Color online) Same as Fig. [I] but J — 2 phase shifts 
and J < 2 mixing parameters are shown. 



In this investigation, we consider three N 3 LO poten- 
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tials which differ by the cutoff parameter A and/or the 
regulator function: 

• A = 414 MeV together with a sharp cutoff (pub- 
lished in Ref. [H]). 

• A = 450 MeV using the regulator function Eq. ([I]) 
with n = 3. We have constructed this potential for 
the present investigation. 

• A = 500 MeV using the regulator function Eq. ([I]) 
with n = 2 for the 2tt exchange contributions. This 
potential was published in 2003 [T]. 

All three potentials use the same (comprehensive) ana- 
lytic expressions which can be found in Ref. 3J. Note 
that the Gaussian regulator function Eq. Q suppresses 
the potential also for Q < A, which is why we use a sharp 
cutoff function in the case of the lowest cutoff of 414 
MeV. Cutoff-independence is an important aspect of an 
EFT. In lower partial waves, the cutoff dependence of the 
NN phase shifts is counter balanced by an appropriate 
adjustment of the contact terms which, at N 3 LO, con- 
tribute in S, P, and D waves. The extent to which cut- 
off independence can be achieved in lower partial waves 
is demonstrated in Figs. [T][2] In F and higher partial 
waves (where there are no NN contact terms) the LECs 
of the dimension-two ttN Lagrangian can be used to ob- 
tain cutoff independence of the phase shift predictions, 
see Table |T] and Fig. [3) 



with qi = pi' — pi, where pi and p/ are the initial and 
final momenta of nucleon i, respectively, and 
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Note that the 2PE 3NF docs not contain any new param- 
eters, because the LECs c±, C3, and C4 appear already in 
the 2PE 2NF. The 1PE contribution is 
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and the 3N contact potential reads 
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In the above, we use gA — 1.29, f v — 92.4 MeV, — 
138.04 MeV, and A x = 700 MeV. The last two 3NF terms 
involve the two new parameters cd and ce, which do not 
appear in the 2N problem. There are many ways to pin 
these two parameters down. The triton binding energy 
and the nd doublet scattering length 2 a n d can be used. 
Alternatively, one may choose the binding energies of 3 H 
and 4 He or an optimal over-all fit of the properties of light 
nuclei. However, in neutron matter, Vd and Ve do not 
contribute such that we do not have to worry about their 
values here. Note also that the C4 term of V c , Eqs.Q and 
([3]), vanishes in neutron matter. 



TABLE I: For the various chiral N 3 LO NN potentials used 
in the present investigation, we show the cutoff parameter A, 
the type of regulator, the exponent n used in the regulator 
function, Eq. dip, and the LECs of the dimension- two tvN 
Lagrangian, c; (In units of GeV -1 ), which are relevant for the 
N 2 LO 3NF in neutron matter. 
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An important advantage of the EFT approach to nu- 
clear forces is that it creates two- and many-body forces 
on an equal footing. The first non-vanishing 3NF occurs 
at N 2 LO. At this order, there are three 3NF topologies: 
the two-pion exchange (2PE), one-pion exchange (1PE), 
and 3N-contact interactions. The 2PE 3N-potential is 
given by 

V = ( 9A Y 1 V • giKgj • gj) pab a b 

C V2/J 2 J^>2 + m2)(^ +TO 2) «* i 

(2) 
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FIG. 3: (Color online) Same as Pig. [TJ but some representative 
peripheral partial waves are shown. 



III. CALCULATION OF THE ENERGY PER 
PARTICLE IN NEUTRON MATTER 

We calculate the ground-state energy per particle 
(g.s.e.) of infinite neutron matter within the framework 
of many-body perturbation theory. In particular, we ex- 
press the g.s.e. as a sum of Goldstone diagrams up to 
third order. 
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FIG. 4: First-, second-, and third-order diagrams of the 
Goldstone expansion included in our calculations with Vnn 
vertices only. Latin-letter subscripts denote particle states, 
greek-letter subscripts correspond to hole states. 
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FIG. 5: Same as in Fig. [4] but also including Vnnn vertices 
which are denoted with thick straight lines. 



In order to take into account the effects of the N 2 LO 
3NF, a density-dependent two-body potential Vnnn is 
added to the chiral N 3 LO potential Vnn- This poten- 
tial Vnnn is obtained by summing one nucleon over 
the filled Fermi sea, which leads to a density-dependent 
two-nucleon interaction [TBI H7j . Hebeler et al. [TU] have 
pointed out that to take care of the correct combinato- 
rial factors of the normal-ordering at the two-body level 
of the 3NF, the matrix elements of V nnnQ^f) are to be 
multiplied by a factor 1/3 in the first-order Hartree-Fock 
(HF) diagram, and by a factor 1/2 in the calculation of 
the single-particle energies (s.p.e.). 

In Fig. [4] we show the diagrams we have included in 
our calculation, where only the Vnn vertices are taken 
into account. The only diagram we do not include is the 
third-order ph diagram. The diagrams that include the 
effects oi Vnnn are shown in Fig. [5] 

The first-order HF contribution is explicitly given by 



The operators P and Q are defined through the rela- 
tionships: 
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In Eq. |7[ the denominator is E(k,k',K) = + 
'"<''•-' •>< ( ' x • I/A. U(k) he- 
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ing the self-consistent single-particle potential: 



The second-order diagrams are computed using the so- 
called angle-average (AA) approximation [TS], and their 
contribution is 
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The particle-particle (pp) and hole-hole (Wi) third- 
order diagrams are also computed in the AA approxi- 
mation, and their explicit expressions are: 
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We have also calculated the [2|1] Pade approximant 



NN potential with a cutoff parameter equal to 500 MeV. 
We have chosen here the potential with the largest cutoff 
since it has the worst perturbative behavior. From the 
inspection of Fig. |6j it can be seen that the energy per 
nucleon calculated at second order, E2, does not differ 
much from the one computed at third order, E 3 , for the 
whole range of densities shown. The perturbative char- 
acter is also indicated by the fact that £3 is quite close 
to the energy obtained with the [2|1] Pade approximant. 
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Ei being the ith order energy contribution in the pertur- 
bative expansion of the g.s.e.. The Pade approximant is 
an estimate of the value to which the perturbative series 
may converge. Thus, the comparison between the third- 
order results and those obtained by means of the [2|1] 
Pade approximant provides an indication of the size of 
the higher-order perturbative terms. 



IV. RESULTS 

As explained in the previous section, we calculate the 
energy per particle of neutron matter in the framework 
of many-body perturbation theory, including contribu- 
tions up to third-order in the interaction. Therefore, it 
is of interest to obtain an idea of the convergence of the 
perturbative expansion of the g.s.e.. 

In Fig. § we show the neutron-matter energy per nu- 
cleon as a function of density, calculated at various orders 
in the perturbative expansion applying the chiral N 3 LO 



FIG. 6: (Color online) Neutron matter energy per particle 
obtained from the N 3 LO 2NF with cutoff A = 500 MeV. The 
first, second, and third order in the perturbative expansion 
and the Pade approximant [2|1] are shown as a function of 
density p. 



We have also investigated the perturbative behavior of 
our calculations when including the effects of Vnnn- Fig. 
[7] shows that, starting from the same N 3 LO potential, 
there is a small enhancement of the higher-order terms 
when including the N 2 LO 3NF. Nevertheless, the results 
at third order are very close to those obtained with the 
[2 1 1] Pade approximant. 

Our main goal is to calculate the g.s.e. per particle 
in infinite neutron matter, starting from N 3 LO chiral 
NN potentials that apply different regulator functions. 
This is done by using the chiral potentials introduced 
in Sec. II. We have added to each 2NF a chiral N 2 LO 
3NF whose low-energy constants ci and C3 are exactly 
the same as in the corresponding N 3 LO iV./V potential, 
see Table H 

In Fig. [8j we show our results, obtained at third-order 
in the perturbative expansion, with and without taking 
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FIG. 7: (Color online) Same as in Fig. [6j but including the 
contribution of the N 2 LO 3NF. 



into account 3NF effects. The results obtained with 2NFs 
show considerable dependence on the choice of the regu- 
lator and its cutoff parameter. This is at variance with 
the desired regulator independence of the EFT. However, 
when including the contributions of the three-body po- 
tentials, which are consistent with their 2NF partner, reg- 
ulator dependence is strongly reduced. This is our main 
result and, at the same time, the first clear evidence that 
modern chiral potentials can provide model-independent 
results in many-body calculations if 2NF and 3NF are 
treated consistently. 

V. CONCLUDING REMARKS AND OUTLOOK. 

In this paper we have studied the regulator dependence 
of many-body predictions when employing chiral two- 
and three- nucleon potentials, using as a testing ground 
the perturbative calculation of the neutron-matter energy 
per particle. We find substantial regulator dependence of 
the predictions when only 2NFs are taken into account. 
The main outcome of this study is the observation that 
the 3NF can play a crucial role in the restoration of reg- 
ulator independence. However, this mechanism works 
properly only when the chiral 2NF and 3NF are treated 
consistently in the sense that the same parameters are 
used for the same vertices that occur in all topologies in- 
volved. This is particularly true for the LECs C\ and C3 
occurring first at N~LO in the chiral power counting. 
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FIG. 8: (Color online) Results obtained for the g.s.e. per par- 
ticle of infinite neutron matter at third-order in perturbation 
theory for three sets of chiral interactions which differ by the 
cutoff parameter A. 



In closing, we note that the present investigation deals 
only with identical nucleon systems, and that the regu- 
lator dependence should also be investigated in systems 
with different concentrations of interacting protons and 
neutrons. In infinite symmetric nuclear matter also con- 
tributions from the intermediate-range l7r-exchange com- 
ponent Vd , and from the short-range contact interaction 
Ve come into play. This means that the calculation of the 
g.s.e. depends also on the coupling constants cd and ce- 
Even though these parameters can be fixed in few-body 
systems, there is some freedom in doing so, resulting in 
more latitude for the 3NF contribution in nuclear matter 
(as compared to pure neutron matter). 

This will be an interesting subject for a future study, 
that may shed more light on the topic of regulator in- 
dependence of many-body calculations with chiral po- 
tentials. The results of such investigations will provide 
valuable guidance for the proper application of these in- 
teractions in microscopic nuclear structure calculations. 
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